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ABSTRACT 



Literature concerning the theory of hit and kill probabilities is 
readily available but is widely scattered. The paper presents a con- 
solidation and a categorization of current literature and provides a 
general discussion of representative models in the field of hit and 
kill probabilities. A literature research matrix is presented to aid 
the researcher in locating existing models which meet his requirements 
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I. INTRODUCTION 



An important problem faced by the weapons system analyist is that 
of determining the probability that a weapons system will effect a 
target hit or kill. The hit or kill probabilities are investigated 
initially by using a single shot or multiple shot hit or kill 
probability models which are further extended to include models which 
consider the entire engagement. Although two factors influence the 
kill probability function - the probability that the weapons system 
can deliver its ordnance to the target and the vunerability of the 
target - only the probability that a given system can deliver its 
ordnance to the target will be discussed in detail in this paper. 

The author f s purpose in presenting a survey of the literature 
reflecting hit and kill probability models is to provide the system 
analyist with a convenient reference for a detailed study using 
models which are pertinent to the weapons system being investigated, 
and where possible to eliminate redundant effort. The majority of 
the models discussed originated in technical memorandum and reports 
from various government agencies or activities and received only 
limited distribution. It is hoped that the consolidation and 
categorization of the models in this paper will enable an analyist 
to more easily locate literature that is relevant to his particular 
needs. Secondly, the paper will serve to familiarize the student or 
researcher with the terminology and assumptions related to hit and 
kill probability models. 

Solution procedures for the models will not be discussed in 
detail as the interested reader can find the solution techniques by 
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reading the original document. The paper focuses upon the logical 
categorization of the models as to the way in which the models are 
formulated and the presentation of relevant models within the 
literature. A convenient categorization of the hit and kill 
probability models is illustrated in figure 1. The literature 
research matrix found in Section III is indicative of the wide range 
and diversification of content of current literature. The terminology 
and symbology varys widely within the literature. To simplify the 
presentation, all models presented in Section IV have been standardized 
so as to conform to the definitions found in Section II. 



6 




Figure 1 
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Figure 1 (continued) 
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II. DEFINITIONS AND TERMINOLOGY 



This section of the paper will present the basic definitions and 
terminology which are used throughout the presentation. Unless other- 
wise specified, the definitions and notations in this section will 
apply to all models and subject material within this paper. 

Hit probability is defined as the probability of a hit or hits 
being made on a target out of a given number of projectiles directed 
at the target. A hit is a blow or impact on a target by a bullet, 
bomb, or other projectile. It should be noted that a projectile 
could technically effect a hit if the force of its explosion created 
a blow or impact to the target. Single shot hit probability is the 
probability that a single projectile fired against a target will hit 
that target under a given set of conditions. Kill probability is 
defined as the conditional probability that a projectile will kill a 
target against which it is fired given that the projectile hits the 
target. As a noun, a kill is that level of damage which destroys a 
target. An aircraft would be considered a kill if the damage was such 
that the aircraft fell out of control immediately after sustaining the 
damage. As a verb, kill means to destroy the target to a specified 
level of damage. Military authorities have designated certain levels 
of damage required for a specific level of kill [Ref. 1]. 

Because of errors inherent in a weapons system the impact point 
for any projectile will in general differ from the intended aim point 
by an amount which is random from projectile to projectile. These 
differences between the intended and actual impact points are the 
result of three basic types of errors. Figure 2 illustrates these 
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Impact Point Error Geometry [Ref. 10] 
Figure 2 
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errors which may be classified into systematic errors, time varying 
errors, and round- to-round dispersion errors. Systematic errors are 
those errors which may be considered constant for the duration of an 
engagement, but can vary from engagement to engagement. Time varying 
errors are errors which vary significantly during the engagement, but 
whose rate of variation is slow relative to the firing rate [Ref. 10]. 

Most models combine systematic and time varying errors into the general 
category of aiming errors. Round- to-round errors are often referred to 
in the literature as ballistic errors which are attributed to the 
ballistics of the projectile and which vary in an uncorrelated manner. 

Bias is a term used when describing errors inherent to a system and 
can result from any of the three error sources. Because bias has an 
explicit meaning within the framework of an individual model, the 
term will be individually defined whenever it is used. 

The distribution of aim points may result from one or a combination 
of factors. Typical factors affecting the distribution of aim points 
are: (a) fire control errors; (b) an unknown target location; (c) 

uncontrolled wander due to the instability of a moving weapons 
carrier; (d) and flight roughness in the case of an aircraft target. 

The aim point may be considered as being fixed or as is the usual 
case, it may be considered as being normally distributed about the 
target center. For the two dimensional case, the aim points are most 
often assumed to have a bivariate normal distribution about the center 
of the target. 

The distribution of rounds about the aim point (ballistic dispersion) 
is random and can be attributed to round imperfections or to meteorological 
conditions. For example, the interaction of temperature on propellents 
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and air density causes uneven force vectors to operate on the pro- 
jectile and thereby causes dispersion of the rounds about the aim 
point. The distribution most often used in describing the round-to- 
round dispersion is the normal distribution. For a large number of 
systems the normal distribution assumption has been verified by 
empirical data. 

Many projectiles inflict damage to a target without having to 
impact the target directly. The blast or explosion effect of a pro- 
jectile or schrapnel from the projectile may cause damage to the 
target if the target is within a lethal distance from the burst point 
of the projectile. A bullet, however, must impact the target to cause 
any degree of damage to the target. Lethality is the term used to 
describe the degree of damage a projectile inflicts on a given target. 
The lethality function is a function which relates projectile-target 
pairs in terms of the distance from the impact point to the projectile. 
The lethality function is the conditional probability of a target kill 
given a specific target is located at point P and the projectile 
impacts at point P-1. The "cookie-cutter” or "zero-one” lethality 
function is used to describe the situation where if the round impacts 
the target or is within a specified distance from the target at impact 
a kill is obtained with probability one; otherwise, the target does 
not sustain any damage. The exponential function is often used to 
model the lethality function of schrapnel producing projectiles. 

The correlation coefficient is a number between plus and minus one 
that indicates the effect of the impact of projectile f i f on the impact 
of projectile ' j 1 , where f i* and 1 j f indicate individual projectiles. 
For independent rounds the correlation coefficient is zero which 
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implies that projectile *i* has no effect upon the impact of projectile 
* j 1 • Similarly, perfect correlation of plus one means that round 1 j 1 
will have the same impact point as round * i 1 • The introduction of the 
correlation coefficient into a mathematical model of hit probability 
is difficult because of the complexity of the resulting function. For 
certain specific models, it can be shown that if an *n' round burst is 
fired, an increase in the correlation coefficient reduces the proba- 
bility of a hit [Ref. 25]. This can be interpreted to mean that if the 
*n' round burst has a fixed aiming error, an increasing number of 
rounds would miss the target as the correlation coefficient was 
increased from zero to one. 

Models are developed which describe single shot and multiple shot 
hit and kill probabilities. The meaning of a single shot event is 
explicit in itself. A multiple shot event can be described in two 
ways. First, such an event is characterized by 1 n* rounds being 
fired with one or more aim points. The rounds may or may not be 
independently aimed; however, the rounds must follow one another in 
sequence (i = l,2,...,n) over a period of time At. A salvo is the 
second type of multiple round event that requires *n* rounds to be 
fired simultaneously at the same target, A salvo of 'n 1 projectiles 
may come from one weapon or may consist of one projectile fired from 
each of 'n 1 weapons. For either event the target must be the same 
for every round of the salvo. The salvo assumption of a single 

V. 

weapon firing 1 n' rounds requires an extremely high rate of fire so 
that the rounds may be assumed to be fired simultaneously. A field 
artillery battery of six guns firing simultaneously at the same target 
exemplifies the salvo fire principle. 
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Problems in the development of models relating hit or kill 
probabilities have two basic conceptual approaches - the coverage 
problem and the vunerable area problem. The coverage problem is 
applicable when the exact location of the target is unknown; but it 
is known that the target is within a specified area. The weapons 
system directed against the target has a lethality function such that 
if the target is within a distance R from the impact point, the target 
will be damaged [Ref. 15]. This method has been used extensively in 
analyzing the effectiveness of a nuclear warhead against a target. 

The vunerable area problem describes the probability of kill in terms 
of damage to one or more vunerable components. The target is reduced 
to one or more vunerable areas (volumes for 3 -dimensional models) 
which if hit, will cause a target kill. The vunerable area is always 
less than the presented target area. If the target has several 
vunerable areas, it is often possible to sum over these areas to get 
an average vunerable area. A typical example of such a target is an 
aircraft which is composed of several vunerable components, one of 
them being the pilot. The aircraft is modeled as a vunerable area 
that is much smaller than the presented area of the aircraft. If 
masking of one area by another occurs, then the problem of target 
representation becomes more difficult. For the two dimensional case, 
an average vunerable area of the target is projected onto a plane 
perpendicular to the trajectory of the oncoming projectile. 

The description of the target in terms of vunerable components 
while being an important element of many models requires a detailed 
analysis of numerous target types with respect to specific types of 
fragments and blast effects. Ballistic Research Laboratories has 
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conducted extensive testing to produce empirical data which forms the 
basis for vunerability input data [Ref. 7]. While the importance of 
such work is not meant to be deemphasized in this paper, the data is 
more relevant to the design of the warhead of a projectile than to 
the weapons system as a whole. The system analyist most often 
utilizes the vunerability data after it has been reduced to a 
vunerable area. 
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III. SUMMARY OF RESEARCH EFFORTS 



To provide a reference guide to the numerous hit and kill 
probability models which exist within the current literature, the 
literature search matrix, figure 3, has been designed. While the 
author does not claim the references listed are a complete listing of 
all available material on the subject, those references listed will 
reduce the set of models the researcher needs to initially investi- 
gate when specific model characteristics are known. Many of the 
references listed include extensive bibliographies which will amplify 
the set of source material available. 

The consolidation and categorization of documented material con- 
cerning subject matter has not been accomplished, or if accomplished, 
is unknown to the author. It was felt that a consolidation of the 
widely scattered material must provide a categorization of the models 
if it were to be of any value. The literature search matrix is the 
vehicle by which the categorization is accomplished. The matrix lists 
six major categories or descriptors of hit and kill probability models. 

By 1 isting the descriptors, a set of applicable models can quickly be 
located by the researcher. The six categories are: (a) dimensionality 

of the model; (b) the number of rounds considered; (c) correlation of 
rounds; (d) lethality function; (e) distribution of aim points; (f) 
and the distribution of the rounds about the aim point. The references 
listed are those found in the bibliography of this paper. To use the 
matrix it is only necessary to locate the characteristic or set of 
characteristics which best describe the problem being investigated. 

The applicable reference or set of references are then easily identified. 
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Th« distribution or rounds about the ala point is assuaed noraal in all nodels except those in Bef 

Tlgure 3 




Throughout the literature there are several basic documents which 



are used continually as source documents for various papers. Reference 
16, COLLECTION OF ARTICLES ON THE THEORY OF FIRING I , appears to be the 
foundation of many of the dispersion optimizing models. The document 
is a translation of a Russian work edited by Kolmogorov in 1948. The 
concepts presented in his document have been studied extensively in this 
country and have been applied when modeling non-nuclear projectiles. 

Two basic documents relating to the coverage problem are Refs. 15 and 
24 which are current through 1963. Reference 15 is an extremely good 
review of the coverage problem as the authors present specific examples 
of coverage problems and the models used in analyzing the problems; 
whereas, Ref. 24 is a bibliography of coverage problems. Grubbs [Ref. 
13] provides basic documentation for offset kill probabilities and 
discusses related works in his paper which was published in 1962. 

These documents are recommended for the individual who is doing 
initial work in the modeling of hit and kill probabilities and is 
desirous of investigating the basic theory for such models. 
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IV. RELATED MODELS 



A. GENERAL 

In this section of the paper models are presented which describe 
the probability of a target hit or kill assuming that a successful 
launch of the projectile has been accomplished. The models are 
categorized by the number of rounds represented within the model - 
single, multiple 'n 1 round, and multiple round salvo fire. Such a 
categorization minimizes the overlap between models and provides the 
researcher with a ready reference for a specific model. If a model 
can readily be extended to another category, it will be so indicated 
The most general model for hit probability is one which is 
equivalent to having no system bias or to having a fixed aim point. 
Let F(r) be the probability density function describing the distri- 
bution of impact points of the projectiles, where r is a symbol for 
the coordinate system representing the location of the impact points 
If 'A 1 denotes the area of the target, then the single shot hit 
probability can be written 

P h = f f(r) dr. (4-1) 

A 

To define the probability of a target kill, a function p(r) is intro 
duced. The function p(r) is the probability that a target is killed 
given the projectile impacts or bursts at r. The probability that a 
projectile kills the target is given by 

P k = J p(r) f (r) dr, (4-2) 

where the integration is taken over the entire burst region [Ref. 7] 
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The distribution of impact points is most often assumed to be 
normal or multivariate normal. For the individual desiring a concise 
review of the bivariate normal distribution, it is recommended that 
Ref. 4 be studied. The mathematical development of the bivariate 
normal distribution was given with special attention being focused 
upon the circular normal distribution. 

B. SINGLE SHOT MODELS 

McNolty has written three papers describing single shot kill 
probability models in which the functions of the weapon lethality, 
weapon bias, and target location are investigated. The lethality 
function of the projectile was modeled in Ref. 18. Expressions for 
the single shot kill probability are derived when the concept of a 
lethal circle of radius R (kill or no kill) cannot adequately repre- 
sent the system being modeled. He defined four functions for the 
weapon lethality (conditional kill probability, P c (r)) whose curves 
are monotonically decreasing from 1 to 0, depending upon the distance 
f r f of the burst from the target. The specific lethality functions 
investigated were 




r > R, 



0 < r < R, 



the zero-one function, 



(4-3) 



P c (r) = exp (-r^/2 b^), (0 < r < ») 



(4-4) 



P c (r) = exp (-b/r), (0 < r < <», b > 0) 



(4-5) 




(0 < r < b , b > 0) 



(r > b) 



(4-6) 




1 - r/b, (0 < r < b, b > 0) 



0. (r > b) 



(4-7) 
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McNolty presented derivations of formulas for the probability of 
killing a point target and for the expected coverage of a population 
of targets when the bias of the weapon system was randomly distributed 
in accordance with a prescribed density function [Ref. 17]. He defined 
weapon bias as the offset distance from the aim point of the center of 
gravity of the weapon impact points. McNolty considered six cases in 
two and three dimensions in which the bias distributions were gamma, 
beta, Maxwe 11-Boltzmann , and Raleigh. Such models are useful if a 
sampling of the weapons system indicates the bias is randomly 
distributed . 

The density of the target was the third basic variable studied by 
McNolty [Ref. 20], He derived expressions for the probability of 
killing a randomly located point target and for the expected coverage 
of an area target of variable density. In this paper McNolty combined 
his two previous works concerning randomly distributed bias and lethal 
effects with the concept of a non-uniform target density. Two specific 
examples were provided which give insight into the use of his models. 

Grubbs provided a useful analytic procedure for computing the 
circular and noncircular offset probabilities of hitting [Ref. 13]. 

He approximated the probabilities of hitting by the use of a central 
chi-square distribution with a fractional number of degrees of freedom 
or a transformation to approximate normality. The paper is signifi- 
cant in the field of offset probabilities of hitting in that he 
discussed related works and also provided six example problems using 
the techniques he presented. His work can be extended to kill 
probability models by incorporating lethality or vunerability data. 

A review of the literature concerning a class of coverage problems 
was prepared by Guenther and Terrangrio in 1963 [Ref. 15]. The problem 
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of killing a point target by a weapons system with a designated 
killing radius was discussed in detail. The paper was divided into 
three sections. The first section was devoted to probability content 
problems in which a spherical region C and a point, B = (b^ ,b^ , . . .b^) , 
would be captured by a sphere of radius R whenever r X r fell within or 
on 

E i=i (x i ' b i )2 = r2 - (4_8) 

Models using the zero one damage function where the target was not a 
point target were discussed in section two of the paper. Section 
three was an extension of sections 1 and 2 where damage functions 
other than the zero one function were considered. This paper not only 
provided modeling techniques, but it also had an extensive biblio- 
graphy which would be useful to the researcher. 

C. MULTIPLE N ROUND MODELS 

When the weapon system analyist must investigate the hit or kill 
probability associated with multiple rounds, mathematical simplicity 
is not usually ameanible to the problem. Except for very specific 
values for the model variables, the mathematical expressions for the 
probabilities often require digital treatment. McNolty derived integral 
expressions for the probabilities for multiple rounds when the target 
was randomly located according to an offset circular normal distri- 
bution and remained in its unknown positions throughout *n* independent 
tosses (shots) of a lethal circle. In his paper, McNolty presented 
integral expressions for the probability of: (1) killing the target 
at least once in r n f tosses (shots) of the lethal circle; (2) killing 
the target exactly T n f times in r N r tosses; (3) requiring less than or 
equal to f m f shots to kill the target exactly once; (4) killing the 
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target at least once in ‘n* tosses when the bias (offset distance) was 
randomly distributed. It is interesting to note that McNolty derived 
an integral expression for killing the target at least once in *n r 
tosses (shots) P^, which was not equal to 

1 - (1 - P R ) n ] (4-9) 

where was the single shot kill probability. His formulation was 
a lways 

< 1 - (1 - P k ) n , (4-10) 



with equality holding only when the target remained in its unknown 
(fixed) location throughout the 'n* shots of the lethal circle. This 
required integration over the entire X,Y plane to account for all 
possible target locations. McNolty also presented the distribution 
for the offset distance , r* (system bias) 



g(r)dr = [2/r(X)]a/a) X r 2X_1 



(4-11) 



•exp (- X/ar ) dr, 

where E(r) = *fa /X * [T (X+l/2)/r (X) ] , and E(r 2 ) = a, 

which included the Raleigh, Maxwell-Boltzmann, and one sided Gaussian 
distributions as special cases. (In most cases for random bias 
application X falls in the range 1/2 < X < 3 [Ref. 19].) 

The diffused target concept was presented by Thompson for multiple 
shot kill probabilities [Ref. 27]. The technique is applicable as an 
event becomes increasingly likely with decreasing miss distance. The 
basic assumption for the concept was that the radial normal distri- 
bution was valid for the miss distance and the weapon bias, if it 
existed, was fixed. A more complete development of the concept was 
found in Ref. 1, where the problem was graphically presented and the 



23 



detailed mathematical model was developed Figure 4 shows the most 
general two dimensional model for the diffused target concept. It 
can be shown that for systems having a fixed bias and a fixed target 
location, the single shot kill probability can be increased by 
increasing the bias up to a specific value. After the bias exceeds 
the lethal radius of the target, however, the kill probability drops 
off rapidly [Refs 1, 25]. 

Breaux investigated the diffused target concept and showed that 
the series solution in one case reduced to the incomplete beta 
function and for a second problem, he derived a new series solution 
[Ref. 9]. Breaux included in his presentation a discussion of com- 
putational considerations as the number of shots increased. It was 
found that the incomplete beta function could be used if the number 
of rounds fired was less than 50; whereas, the series Breaux developed 
was valid for up to 1,000 rounds. 

The extension of a multiple shot model into an engagement kill 
model was done by Scheu in a paper investigating aircraft attrition 
by gun systems [Ref. 26]. In his development, he assumed that the 
single shot kill probability was a function of slant range, and then 
he averaged the single shot kill probability over the aircraft pass 
through the killing zone of the weapon. The engagement kill proba- 
bility was developed to include the effect of terrain masking. The 
mathematical models were validated with empirical data. The validation 
of the model was not surprising since the empirical data was initially 
used in the assignment of specific values to constant factors which 
were part of the mathematical model. The target was represented by 
an average vunerable area and the engagement kill probability was 
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(Pr) is single-shot kill probability 

SSK 

h, k are bias components 

a is vunerable radius 
c 

is standard deviation of burst 
pattern dispersion 

2 

a is vunerable area = 2 tt a 

c 



Two Dimensional Fixed Bias Burst Pattern 
Figure 4 
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EXP = 1 - exp(-Mp), (4-12) 

where Mp was a function of the number of rounds fired, terrain, air- 
craft velocity, single shot kill probability, and the slant range 
between the gun and the aircraft. 

The models presented thus far have all assumed that the rounds 
were independent. Peterson developed a method for finding the 
conditional hit and miss probabilities on a rectangular target given 
a hit or a miss on the previous shot when part of the random error was 
common to both shots but varied from engagement to engagement. The 
correlation coefficient depended upon the relative size of the two 
elements of the total dispersion. He provided two methods of achieving 
computational as well as tabular results [Ref. 6]. 

The effect of correlation between rounds was modeled by Helgert 
where consideration was also made for time varying errors [Ref. 10]. 
Although a detailed development of the model was not presented, the 
model was used in a detailed example with empirical data. Graphs and 
tables of cumulative hit probabilities were provided which showed the 
effect of a change in the error statistics for a given number of 
rounds fired (10 to 200). The target was assumed to be nearly 
circular or square, no aim point adjustment during the engagement, and 
equal error statistics in both coordinates. The error statistics 
considered were - target area, number of rounds fired, standard 
deviation of round-to-round dispersion errors, standard deviation of 
systematic errors, standard deviation of time varying errors, and the 
"average correlation between rounds". 
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D. MULTIPLE ROUND SALVO MODELS 



Weapons systems which deliver a salvo of rounds to a target require 
a unique solution process to solve the kill probability problem. The 
most general assumption regarding the salvo model is to assume that 
the salvo aim point is distributed around the target according to a 
bivariate normal distribution with the individual rounds being subject 
to the bivariate normal distribution about the center of impact. Using 
these two basic assumptions, Rice and Bottero developed detailed kill 
probability formulae which can be conveniently used for hand compu- 
tation [Ref. 22]. The weapon lethality was specified through the use 
of a near miss zone surrounding the two and three dimensional target 
and a conditional kill probability. The most general model in the work 
included the target orientation as a variable. Because of the clarity 
of the material and the logical presentation of the formula derivation, 
the model will be discussed in detail as an excellent example of the 
approach to salvo kill probability problems. The general derivation 
was : 

a. Develop an expression for the probability that a single round 
of the salvo will impact within the target [p(x,y)]. 

b. Using P^^*p(x,y) as the probability that the round will kill 

the target where P ^ is the conditional probability of a kill 
given a target hit, the probability that the target will 
survive the round is 1 - *hk P r °h a hility that 

the target survives all , n l rounds of the salvo is then 

[1 - P hk -p(x,y)] n . 

c. The probability of success A(x,y) can then be written as 
1 - [1 - *hk which can be interpreted as the 
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probability that at least one round of the salvo with aim 
point (x,y) will produce the desired level of damage. 

d. The probability P^(n) that at least one damaging hit will be 
made by a salvo of , n l rounds with intended aim point (x^y 1 ) 
is derived by averaging A(x,y) over all possible x and y. 

e. The probability P^(n) °f at least q damaging hits in a salvo 
of , n l rounds can also be developed. 

A computer program has been written for three dimensional targets 
which had the orientation angle assumed to be zero. The integrals 
were evaluated by using a Hasting approximation formula and Wedde 1 s 
rule . 

A comparative study of three mathematical models used in esti- 
mating salvo hit probabilities was conducted by Kline and Fergerson 
[Ref. 23]. The first model assumed that the landing position of the 
individual rounds of the salvo were independently and identically 
distributed. The second model assumed that the salvo of weapons was 
dispersed in a well defined rectangular array. The third model 
investigated assumed that the weapons which formed the salvo were 
randomly distributed within a randomly distributed circular region 
where the center of the circle was determined by a random selection 
from a circular normal distribution. The significant difference 
between the three models was the assumed statistical dependence 
between the landing positions of the various weapons of the salvo. 

The landing position of each weapon was assumed to be a circular 
normal distribution function in each of the three models. Strict 
statistical dependence between weapons was exemplified by model two 
with strict independence evident in model one. Model three repre- 
sented a mixture of the two extremes. The paper graphically depicted 
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the fact that the hit probabilities determined were directly affected 
by the dependence criteria (model) selected. 

Banish developed a model for salvo kill probability using matrix 
algebra [Ref. 8]. His model took into consideration the time 
variational error within the system. His development assumed a 
normal distribution of the aim point and a normal distribution of 
the rounds about the aim position with a correlation between rounds 
being present. Banish r s study was unique in that he utilized matrix 
algebra in the development of his model. 

An expansion of the work done by Kolmogorov [Ref. 16] was accom- 
plished by Walsh in his study of salvo kill probability models [Refs. 
28, 29]. Walsh developed a general expression of the salvo kill 
probabilities which could be approximated using an explicit function 
of four parameters. The four parameters were the number of rounds 
'n f , the average target vunerability times the target size ’VT 1 , the 
salvo aiming error dispersion f a^ f , and a function of the round 

dispersion 1 q 1 and two dispersion type terms which depended on the 
R 

2 2 

target size, shape, and vunerability 1 G,^, • The probability 

that exactly 'm' rounds of the salvo resulted in target kills could 
therefore be expressed as 




rC-l) 1 (VT) wri ' 1 n!l/fi!m!(n-m-i). , (2n) : 



n!]/[i!rn!(n-m-i)!(2n) 1 ^ i ] 




(4-13) 



The terms G,^ and were dependent upon the target size, shape, and 



vunerability such that 




(4-14) 



T 
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(4-15) 



where V(u,v) was the conditional probability of a target kill given 
the round hits the target at (u,v) and V was the average target 
vunerability [Ref. 28]. Walsh also presented a procedure for evalu- 
ating the maximum salvo kill probability obtainable and for deter- 
mining the optimum round hit- location-probabi lity distribution that 
yields that maximum value. The method of approach to the problem of 
maximizing the salvo kill probability required the application of the 
calculus of variations. The determination of the approximately 
optimum form of the probability density function necessitated an 
iterative *cut-and-try 1 method. The maximum value of the salvo kill 
probability P^ was 



where A(x) was the density function for the probability distribution 
of the aiming error x and p(x) was the optimal form of the con- 
ditional probability that a round with expected hit location (x) had 
enough individual effect to kill the target. The method developed 



where V was the average target vunerability and T was the size of the 
target (area for d = 2, volume for d = 3). The integration was over 
values of x such that A(x) > C. The maximum salvo kill probability 
was expressed as 




A(x)dx. • • -dx , (4- 16) 

i d 



for evaluating P^ and p(x), (n > 1), was to first determine the value 
of C from 




(4-17) 
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= 1-C ^[size of x-region for which A(x) > C] - VT 

/••/ A (x)dx ^ • • -dx^, (4-18) 

A(x) < C 

f 1- [C/A(x) ] ^ ^° ^ x such that A(x) > C 
p(x) = < (4-19) 

0 otherwise. 

Walsh presented a numerical example to illustrate the 1 cut-and-try 1 
method for determining the optimum form for the probability density 
function of the round hit locations [Ref. 29]. 

A computational procedure for estimating the salvo hit proba- 
bilities for offset circular targets was developed by Groves and Smith 
in 1956 [Ref. 12]. The approach taken by the authors was to express 
the distribution of the impact (burst point) of the missiles in terms 
of the distance from the center of the target rather than from the 
center of impact of the salvo when the center of impact of the 
missiles was offset from the center of the target. Sample problems 
were illustrative of the proposed methodology and graphs were utilized 
to aid in the solution procedure. 

Grubbs expanded his earlier work on circular and noncircular 
offset probabilities to include the development of an expression for 
the expected fraction of damage to a target [Ref. 14]. In general 
the average fraction of damage over a circular target was 




7<n> - // 



h(u,v) J(u,v) du dv, 



2 

u 



+ V 






(4-20) 



where h(u,v) was the distribution of the target elements over a 
circular target of radius R T ^ whose center was at the origin and 
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J(u,v) was the probability that at least one of the 'n 1 rounds of the 
salvo would damage a target element at (u,v). The final result after 
several approximations was 



f(n) - [ 2 ^'V + ^)(o? • ft, y )/R ?] Cl CD <- 1)1+1 - 1/1 



J 2 2 , f ( 2 . 2. 2 2 V 

• <a, a, A/ (a, +a )(a. +a ), 

: kx ky kx x ky y " 

1-exp [-i R^/2 (a^+4+iaJ)]l 



1/2 



• < 1-exp [-i R^/2 (a^ y +ay+iaj)][ 



1 

i] 1/2 






(4-21) 



where a, and a, were the standard deviations of the noncircular 
kx ky 

damage function. The above formula was a general formula for the 
noncircular case and uniform target density. A computational pro- 
cedure for the binomial coefficients using Jacobi Polynomials was 
developed by Breaux and Mohler which simplified the evaluation for 
the approximation [Ref. 3]. 

The development of an approximate formula for salvo kill proba- 
bilities by Walsh [Ref. 28] can be contrasted to the expected 
fractional damage approximation by Grubbs [Ref. 14]. Walsh used a 
general lethality function and the concept of target variances com- 
bined with the idea that the round to round variability strongly 
dominated the target variances. Circular normal delivery errors 
were also a basic assumption. Grubbs, however, used a lethality 
function that was treated as an elliptical normal fall-off type 
pattern. Noncircular delivery errors could easily be included in the 
theory developed by Grubbs. 
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For the analyist not having access to a high speed computer. 
Groves developed a method for the hand computation of the expected 
fractional kill f F f for an elliptical target by a number of rounds 
! n' all of which have the same lethal area 'A 1 and are all aimed at 
the same point (x,y) relative to the target center [Ref. 5]. By 
letting a and b represent the semi-axis of the target with and 
being the delivery error standard deviations, the expected fractional 
kill could be denoted 



i -1 



F(a,b,x,y,a x ,C y) A,n) = — 



where 



CCDG&) 

r~A ■ 2 „ r~A i 2~ , 

a "/s + V 13 l2^ + V ( 22) 



is a function that gives the probability of hitting a circular 

target of radius a/I centered at the origin of a coordinate system 

given a bivariate normal distribution of hits with mean and 

a 

standard deviation ~ in one direction: and mean and standard 

a b _ 

deviation ^ in the other direction. The values for P. f*/i ,— , 7 * 

b hV a b a 5 b y 

must be determined using approximation formulae or tabular values. 
Groves referenced both formulae and tables which may be used by the 
analyist and included tabular values for for the special case when 
x = y = 0 . 



Because of the difficulty involved in using the series developed 
by Groves for the expected fractional kill of a circular target when 
the number of rounds 'n f is large, Breaux proposed an alternate series 



f(n) = C [(1-xe‘V - (l-x) 1 ]/!. (4-23) 

The terms of the series decrease monotonically with the summation 
index i and were therefore readily summable [Ref. 2], 
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The models presented in this section are representative of the 
approaches which have been taken in hit and kill probability models. 

The models represent coverage problems, integral approximations, 
diffused target concepts, and techniques of numerical approximations. 
The relevancy of a specific model will be dependent upon the character- 
istics of the weapon system being analyzed. It was therefore felt that 
the detailed examination of each model was unnecessary and that the 
generalized presentation would give an insight into the major 
techniques and assumptions associated with hit and kill probability 
mode Is . 
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V. AREAS FOR FUTURE RESEARCH 



Areas for future research in hit and kill probability modeling are 
open both in the theoretical development of the models and in the 
application of existing models for the evaluation of specific weapons 
systems. Specific recommendations for future study include: 

a. The optimization of the distribution of rounds about the aim 
point for salvo fire. Such a study should include the intro- 
duction of artificial dispersion. 

b. The effect of time varying errors upon the kill probability by 
considering the relative movement between the target and the 
projectile . 

c. The effect of round-to-round correlation upon the hit or kill 
probability for multiple round models. 

d. The dynamic modeling of hit and kill probabilities when the 
probabilities are a function of time or range from the weapon 
to the target. 

e. The development of a dynamic model that reflects the effect of 
varying the rate of fire of the weapons system upon hit or kill 
probabilities . 
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